Substrate Specificity of Peptide Adsorption: A Model Study 
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Applying the contact density chain-growth algorithm to lattice heteropolymers, we identify the 
conformational transitions of a nongrafted hydrophobic-polar heteropolymer with 103 residues in 
the vicinity of a polar, a hydrophobic, and a uniformly attractive substrate. Introducing only two 
system parameters, the numbers of surface contacts and intrinsic hydrophobic contacts, respectively, 
we obtain surprisingly complex temperature and solvent dependent, substrate-specific pseudo-phase 
diagrams. 
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From recent experiments of the adsorption of short 
peptides at semiconductor substrates it is known that dif- 
ferent surface properties (materials such as Si or GaAs, 
crystal orientation, etc.) as well as different amino acid 
sequences strongly influence the binding properties of 
these peptides at the substrate [H, Q. This specificity 
will be of particular importance for future sensory de- 
vices and pattern recognition Q at the nanometer scale. 
The reasons for this binding specificity are far from being 
clear, and it is a big challenge from the experimental and 
theoretical point of view to understand the basic prin- 
ciples of substrate-peptide cooperativity. This problem 
can be seen as embedded into a class of similar studies, 
where the adsorption and docking behavior of polymers 
is essential, e.g., protein-ligand binding prewetting 
and layering transitions in polymer solutions as well as 
dewetting of polymer films molecular pattern, elec- 
trophoretic polymer deposition and growth 

The experimental equipment has reached such a high 
resolution allowing for precise identification of single 
molecule shapes at the substrate, and the available com- 
putational capacities combined with sophisticated algo- 
rithms will make it possible to examine the problem of a 
hybrid interface between biological and inorganic mate- 
rials [7| step by step. 

In a first step, we have recently analyzed the adsorp- 
tion of a finite, nongrafted homopolymer at an attractive 
substrate in a cavity and discussed in detail the phase 
diagram in the thermodynamic limit, as well as pseudo- 
transitions that depend strongly on the given number of 
monomers Similar studies of grafted polymers are 
reported in Refs. [1, [13, HI- For heteropolymers the 
thermodynamic limit is unreachable because of the se- 
quence of different types of monomers. This "disorder- 
inducing" sequence renders the heteropolymer adsorp- 
tion a distinguishingly different problem to homopolymer 
substrate-binding. The pseudo-phase transitions of the 
heteropolymer system will strongly depend on three main 
points: the sequence, the monomer-specific interaction 
with the substrate, and the total number of monomers 
in the chain. In this work, we particularly focus on the 
substrate-specificity. In order to reduce the complexity 



of the problem to a minimum, we study the hydrophobic- 
polar (HP) model [l^, where the heteropolymer consists 
of a given sequence of only two types of monomers: hy- 
drophobic (H) and polar (P). We use the simplest form 
of the model, where only the hydrophobic force acts 
and the number of nearest-neighbor contacts between H 
monomers being nonadjacent along the chain, riHH, is 
related to the energy of the heteropolymer. The interac- 
tion with the substrate is modeled in a like manner: The 
energy of the heteropolymer is reduced by the number 
of nearest-neighbor contacts between the substrate and 
those monomers that experience the attractive force of 
the substrate. For all other monomers the influence of 
the substrate is only entropic. 

In order to study the specificity of surface-binding, we 
investigate three attractive substrate models. In the first 
variant, all monomers, independent of their hydrophobic 
or polar character, are equally attracted by the substrate 
and the energy of the system is proportional to the to- 
tal number of monomer-surface contacts, ri^^^. In the 
second and third model, the substrate is either hydropho- 
bic or polar, i.e., in the first case only the hydrophobic 
monomers in the heteropolymer sequence are attracted, 
and in the latter the attraction between substrate and 
heteropolymer dipoles dominates. In these models, the 
respective hydrophobic substrate contacts, n^, and polar 
surface contacts, , are energetically favored. Thus, the 
models can be expressed in the energetic form 

Es{ns,nuii) = -esUs - eKHnnH, (1) 



where, depending on the substrate model, Ug = ii^^^ , 
n^, or . For our qualitative study, it is sufficient 
to choose for all three models the same energy scales 
£3 = 1 and £hh = s, where s denotes the solubility which 
controls the solvent quality (the larger the value of s, 
the worse the solvent). Since we are interested in the 
fluctuations of the respective contact numbers with re- 
spect to temperature T and solubility s, we define the 
contact density as 5(7is,nHH) = Sn,o g^inim) + (1 — 
Sn^o)g^{ns,nYiM), with the contributions of the densities 
of conformations without (g") and with (g^) contact to 
the substrate. In order to regularize the influence of the 
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unbound conformations and for computational efficiency, 
the heteropolymer is restricted to reside in a cage, i.e., in 
addition to the physically interesting attractive surface 
there is a steric, neutral wall parallel to it in a distance 
Zw The value of is chosen sufficiently large to keep 
the influence on the unbound heteropolymer small (in 
this work we used = 200). Introducing the parti- 
tion sum Zt,s ~ S„,,„HH ^("■si "•hh) ex.p{-Es/kBT) and 
denoting thermodynamic expectation values of a quan- 
tity 0(ns,nHH) by (O), the contact correlation matrix 
Mxy{T,s) = {xy)c = (xy) - {x){y) with x,y = nHH,"-s 
separates the fluctuations of the surface and hydrophobic 
contacts according to the respective energy scale vector 
(ehHj^s) = (s,!), and therefore the energetic fluctua- 
tions are accounted for in the specific heat, defined by 
(from now on we set fc^ = 1) 



Cv{T,s)^^isA)M{T,s) 



(2) 



The heteropolymer sequence we chose in our study pos- 
sesses 103 monomers (37 being hydrophobic, 66 polar) 
as introduced in Ref. [13 [and often used for benchmark 
tests of new algorithms 13, Q, 0, ■ The advantage 
is that this heteropolymer forms a nicely compact hy- 
drophobic core in the low-energy conformations (as de- 
picted for s = 1 and Ug = 0, e.g., in Ref. completely 
screened from the solvent by a shell of polar monomers. 

In order to calculate the contact densities for the three 
systems, we have applied an enhanced version of the mul- 
ticanonical chain-growth algorithm 17, T^. In contrast 



to move-set based Metropolis Monte Carlo or conven- 
tional chain-growth methods which would require many 
separate simulations to obtain results for different param- 
eter pairs (T, s) and which frequently suffer from slowing 
down in the low-temperature sector, our method allows 
the computation of the complete contact density for each 
system within a single simulation run. Since the contact 
density is independent of temperature and solubility, en- 
ergetic quantites such as the specific heat can easily be 
calculated for all values of T and s (nonenergetic quanti- 
ties require accumulated densities to be measured within 
the simulation, but this is also no problem). 

In Figs. [IKa)-(c) the contour profiles of the specific 
heats for the different substrates are shown (the brighter 
the color the larger the value of Cy). We interpret the 
ridges (for accentuation marked by white and gray lines) 
as the boundaries of the pseudo-phases. It should be 
noted, however, that in such a finite system the exact po- 
sitions of active regions exhibited by fluctuations of other 
quantities usually deviate, but the qualitative behavior is 
similar [lij . The gray lines indicate the main transition 
lines, while the white lines separate pseudo-phases that 
strongly depend on speciflc properties of the heteropoly- 
mer, such as its exact number and sequence of hydropho- 
bic and polar monomers. As a flrst result, we have found 
that the binding-unbinding transition appears to be first- 
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FIG. 1: (Color online) Pseudo-phase diagrams of the lOSmer 
near three different substrates that are attractive for (a) all, 
(b) only hydrophobic, and (c) only polar monomers. 



order like. Assuming the contact numbers Hs and nun 
to be kind of order parameters adequately describing the 
state of the heteropolymer, we define the contact free en- 
ergy as Jt,s("s,"-hh) = -Tin [g(ns, tihh) exp(-£;s/r)] 
and the probability for a macro-state with Us sub- 
strate and nHH hydrophobic contacts as pt,s ("-s, ?^hh) = 
g{ns,n}jii)Gxp{—Es/T)/ZT,s- Close to the binding- 
unbinding transition, adsorbed and desorbed states co- 
exist. This is exhibited by two clearly separated minima 
of the contact free energy Ft,s('^s, ^^hh)- In the figures 
we have marked the coexistence line, where both min- 
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ima take the same value, by the dashed black lines. At 
lower temperatures, the most probable conformation is 
an adsorbed one, while for higher temperatures desorbed 
conformations dominate. 

Despite the surprisingly rich and complex phase be- 
havior there are main "phases" that can be distinguished 
in all three systems. These are separated in Figs. [Ha)~ 
(c) by gray lines. Comparing the three systems we find 
that they all possess pseudo-phases, where adsorbed com- 
pact (AC), adsorbed expanded (AE), desorbed compact 
(DC), and desorbed expanded (DE) conformations dom- 
inate, similar to the generic phase diagram of the ho- 
mopolymer 0, ^ [l^, 11 1 . "Compact" here means that 
the heteropolymer has formed a dense hydrophobic core, 
while expanded conformations exhibit dissolved, random- 
coil like structures. The sequence and substrate speci- 
ficity of heteropolymers generates, of course, a rich set of 
new interesting and selective phenomena not available 
for homopolymers. One example is the pseudo-phase 
of adsorbed globules (AG), which is noticeably present 
only in those systems, where all monomers are equally 
attractive to the substrate (Fig. [IJa)) and where polar 
monomers favor contact with the surface (Fig. [Tl^b)). In 
this phase, the conformations are intermediates in the 
binding-unbinding region. This means that monomers 
currently desorbed from the substrate have not yet found 
their position within a compact conformation. There- 
fore, the hydrophobic core, which is smaller than in the 
respective adsorbed phase (i.e., at constant solubility s), 
appears as a loose cluster of hydrophobic monomers. 

In Figs. [2l^a)-(c), we have plotted, exemplified for 
s — 2, the statistical averages of the contact numbers 
Us and nHH as well as their self- and cross-correlations 
M for the three systems. For comparison we have also in- 
cluded the specific heat, whose peaks correspond to the 
intersected transition lines of Figs, ma)-(c) at s = 2. 
From Figs. Elja) and (c) we read off that the transition 
from AC to AG near T « 0.4 is mediated by fluctuations 
of the intrinsic hydrophobic contacts. The very dense hy- 
drophobic domains in the AC subphases lose their com- 
pactness. This transition is absent in the hydrophobic- 
substrate system (Fig. [2jb)). The signal seen belongs to 
a hydrophobic layering AC subphase transition, which in- 
fluences mainly the number of surface contacts n^. The 
second peak of the specific heats belongs to the tran- 
sition between adsorbed compact or globular (AC, AG) 
and expanded (AE) conformations. This behavior is sim- 
ilar in all three systems. Remarkably, it is accompanied 
by a strong anti-correlation between surface and intrinsic 
contact numbers, Us and riHH- Not surprisingly, the hy- 
drophobic contact number nmi fluctuates stronger than 
the number of surface contacts, but apparently in a dif- 
ferent way. Dense conformations with hydrophobic core 
(and therefore many hydrophobic contacts) possess a rel- 
atively small number of surface contacts. Vice versa, con- 
formations with many surface contacts cannot form com- 
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FIG. 2: (Color online) Temperature dependence of specific 
heat, correlation matrix components, and contact number ex- 
pectation values of the lOSmer for surfaces attractive for (a) 
all, (b) only hydrophobic, and (c) only polar monomers at 
s = 2. 



pact hydrophobic domains. Finally, the third specific- 
heat peak marks the binding-unbinding transition, which 
is, as expected, due to a strong fluctuation of the surface 
contact number. 

The strongest difference between the three systems 
is their behavior in pseudo-phase AC, which is roughly 
parameterized by s > 5T. If hydrophobic and po- 
lar monomers are equally attracted by the substrate 
(Fig.[lja)), we find three AC subphases in the parameter 
space plotted. In subphase ACl, film- like conformations 
dominate, i.e., all 103 monomers are in contact with the 
substrate. Due to the good solvent quality in this region, 
the formation of a hydrophobic core is less attractive than 
the maximal deposition of all monomers at the surface, 
the ground state being (n^"*"^, nHH)min = (103,32). In 
fact, instead of a single compact hydrophobic core there 
are nonconnected hydrophobic clusters. At least on the 
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used simple cubic lattice and the chosen sequence, the 
formation of a single hydrophobic core is necessarily ac- 
companied by an unbinding of certain polar monomers 
and, in consequence, an extension of the conformation 
into the third spatial dimension. In fact, this happens 
when entering AC2 [(nf nHH)min = (64,47)], where a 
single hydrophobic two-layer domain has formed at the 
expense of losing surface contacts. In ACS, the het- 
eropolymer has maximized the number of hydrophobic 
contacts and only local arrangements of monomers on 
the surface of the very compact structure lead to the still 
possible maximum number of substrate contacts. -Pr.s is 
minimal for (n^"'"^, nHH)min = (40,52). 

The behavior of the heteropolymer adsorbed at a sur- 
face that is only attractive to hydrophobic monomers 
(Fig.[ljb)) is apparently different in the AC phase. Since 
surface contacts of polar monomers are energetically not 
favored, the subphase structure is determined by the con- 
currence of two hydrophobic forces: substrate attraction 
and formation of intrinsic contacts. In ACl, the num- 
ber of hydrophobic substrate contacts is maximal for 
the single hydrophobic layer, (n™,?iHH)min — (37,42). 
The single two-dimensional hydrophobic domain is also 
maximally compact, at the expense of displacing polar 
monomers into a second layer. In subphase AC 2, intrinsic 
contacts are entropically broken with minimal free energy 
for 35 < nHH < 40, while — 37 remains maximal. 
Another AC subphase, AC3, exhibits a hydrophobic lay- 
ering transition at the expense of hydrophobic substrate 
contacts. Much more interesting is the subphase tran- 
sition from ACl to AC5. The number of hydrophobic 
substrate contacts n™ of the ground-state conforma- 
tion dramatically decreases (from 37 to 4) and the hy- 
drophobic monomers collapse in a one-step process from 
the compact two-dimensional domain to the maximally 
compact three-dimensional hydrophobic core. The con- 
formations are mushroom-like structures grafted at the 
substrate. AC4 is similar to ACS, with advancing des- 
orption. 

Not less exciting is the subphase structure of the het- 
eropolymer interacting with a polar substrate (Fig.[T]Jc)). 
For small values of s and T, the behavior of the het- 
eropolymer is dominated by the concurrence between po- 
lar monomers contacting the substrate and hydrophobic 
monomers favoring the formation of a hydrophobic core, 
which, however, also requires cooperativity of the polar 
monomers. In ACl, film-like conformations (n^ = 66, 
riHH = 31) with disconnected hydrophobic clusters dom- 
inate. Entering AC2, hydrophobic contacts are energet- 
ically favored and a second hydrophobic layer forms at 
the expense of a reduction of polar substrate contacts 
[(n^, nHH)min — (61,37)]. In AC3, the upper layer is 
mainly hydrophobic [(nf , nHH)min = (53,45)], while the 
poor quality of the solvent (s large) and the compara- 
tively strong hydrophobic force let the conformation fur- 
ther collapse [AC4: (nf,nHH)min = (42,52)] and the 



steric cooperativity forces more polar monomers to break 
the contact to the surface and to form a shell surrounding 
the hydrophobic core [AC5: (n^ ,nHH)min — (33,54)]. 

Summarizing, we have performed a detailed analysis of 
the pseudo-phase diagrams in the T-s plane for a selected 
heteropolymer with 103 monomers in cavities with an 
adsorbing substrate being either attractive independently 
of the monomer type, or selective to hydrophobic or polar 
monomers, respectively. Beside the expected adsorbed 
and desorbed phases, we find a rich subphase structure in 
the adsorbed phases with compact hydrophobic domains, 
which is specific to heteropolymers. In particular, the 
formation of layered subphases in the low-temperature 
region depends mainly on the quality of the solvent. 

Here, we have mainly focused on the contact numbers 
ris and nnH, but the study of structural quantities, such 
as the gyration tensor, which exhibits a phase-dependent 
asymmetry in the components parallel and perpendicu- 
lar to the substrate, confirms our interpretation of the 
subphase-behavior of the system [l^ . Since current 
experimental equipment is capable to reveal molecular 
structures at the nanometer scale, it should be possible 
to investigate the grafted structures dependent on the 
solvent quality. This is essential for answering the ques- 
tion under what circumstances binding forces are strong 
enough to refold peptides or proteins. The vision of fu- 
ture biotechnological and medical applications is fasci- 
nating as it ranges from protein-specific sensory devices 
to molecular electronic devices at the nanoscale. 
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